Complete QED theory of multiphoton trident pair production in strong laser fields 
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Electron-positron pair creation by multiphoton absorption in the collision of a relativistic electron 
with a strong laser beam is calculated within laser-dressed quantum electrodynamics. Total produc- 
tion rates, positron spectra, and relative contributions of different reaction channels are obtained in 
various interaction regimes. We study the process in a manifestly nonperturbative domain which is 
shown accessible to future experiments utilizing the electron beamlines at novel x-ray laser facilities 
or all-optical setups based on laser acceleration. Our theory moreover allows us to add further 
insights into the experimental data from SLAC [D.Burke et al, Phys. Rev. Lett. 79, 1626 (1997)]. 
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Quantum electrodynamics (QED) describes the inter- 
action of electrons and photons. When the photon source 
is an intense laser field, electrons may couple nonlinearly 
to the field giving rise to multiphoton processes. Theo- 
reticians started to consider nonlinear QED phenomena 
such as multiphoton Compton scattering or e + e~_pair 
production soon after the invention of the laser [l], [2|. 

In the mid 1990s pioneering studies at SLAC (Stanford, 
California) revealed nonlinear QED effects in experiment. 
In particular, the first and so far unique observation of 
multiphoton pair production was accomplished in colli- 
sions of the fa 50 GeV electron beam from SLAC's linear 
accelerator with an intense laser pulse [3|. The laser fre- 
quency and field strength were largely Doppler-enhanced 
in the rest frame of the high-energy projectile. 

The SLAC experiment has triggered substantial theo- 
retical efforts on laser-induced e + e~ pair creation in the 
last decade (see @,|^] for reviews). Most studies consider 
pair creation in laser-proton collisions (e.g., |q -lll|) or in 
counterpropagating laser beams (e.g., |ia4l7{ ). Mean- 
while, theoreticians are already considering more refined 
aspects of the process such as final-state pai r correlations 

[§] , the influence of more complex fields [111 Hij , and the 
+ , 



creation of fi r \i 



pairs . 

Nevertheless, despite all these efforts, no complete 
QED calculation of the SLAC experiment exists as yet. 
One purpose of the present study is to fill this gap in the 
theory of strong-field phenomena, taking advantage of 
the recent theoretical progress and advancing it further. 
Indeed, nonlinear pair creation, both in the perturbative 
and nonperturbative domains, is nowadays becoming ac- 
cessible to all-optical setups involving laser acceleration 
devices. Note that a recent experiment on powerful laser- 
solid interaction observed abundant (linear) pair creation 
through the conversion of 7-ray bremsstrahlung in the 
field of gold nuclei [jjj . 

In comparison with multiphoton pair creation in laser- 
proton collisions, the theoretical consideration of laser- 
electron collisions is rendered more involved in several 
respects. One needs to take into account (i) the dress- 



ing of the projectile electron by the laser field, (ii) the 
recoil the projectile suffers during the collision, (iii) the 
indistinguishability of the scattered projectile with the 
created electron of the pair (Pauli principle), and (iv) 
the possibility of real photon emission by the projectile 
in the field (Compton scattering). 

In laser-electron collisions, two pair creation processes 
are usually distinguished. The first is of Bethe-Hcitlcr 
type; the pair is produced by the absorption of N laser 
photons in the Coulomb field of the incoming electron: 



e + Nlo -> e' - 



e e 



(1) 



The second is a two-step process where first a high-energy 
7-photon is generated by Compton backscattering off the 
electron beam, which afterwards creates the pair in a 
photon-multiphoton collision [l|, Q : 



7 + No; -> e H 



(2) 



Reaction ([2]) represents the strong-field generalization of 
the process 27 — » e + e~ first studied by Breit and Wheeler 
[l9l | and exhibits a nonperturbative nature at very high 
fields. We name it as multiphoton Breit- Wheeler process. 

The analysis of SLAC's experimental data 0] relied on 
separate simulations of both processes. However, while 
for the multiphoton Breit- Wheeler reaction a sophisti- 
cated nonperturbative theory was available [H, Q], the 
contribution from the direct process was estimated in a 
rather approximate manner [201 ] based on the Weizsacker- 
Williams method since a formal theory did not exist. 

In this Letter we provide a nonperturbative laser- 
dressed QED calculation of multiphoton trident pair cre- 
ation in laser-electron collisions. Our approach treats 
the competing processes ([1} and © in a unified way and 
opens deeper insights into the SLAC measurements. Fur- 
ther, we evaluate the creation rates in the fully nonper- 
turbative regime which could be probed by future exper- 
iments employing upcoming technologies. 

The laser-dressed Feynman diagrams of multiphoton 
trident pair creation are depicted in Fig.[T] The leptons 
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FIG. 1: Furry-Feynman diagrams of multiphoton trident pair 
production in electron-laser collisions. The zigzag-lines repre- 
sent the exact lepton wave-functions in the laser field (Dirac- 
Volkov states [21() and are labeled by the laser-dressed parti- 
cle momenta. In the left diagram, the incoming electron scat- 
ters from a state of dressed momentum q to q by absorbing 
n laser photons and emitting an intermediate photon, which 
afterwards decays into an e + e~ pair upon absorption of n' 
laser photons. The corresponding exchange diagram is shown 
on the right. In the weak-field limit (£ <§C 1), the diagrams 
can be expanded in Feynman diagrams of field- free QED. For 
a total number of N = n + n' — 1 absorbed photons, eight 
leading-order diagrams arise this way [52]. The case N = 6 
corresponds to 168 such diagrams. 



are described by Dirac-Volkov states 21] which fully ac- 
count for their interaction with the external plane-wave 
laser field (Furry picture). To leading order in the QED 
finestructure constant a the amplitude for the process 
reads Sfi = M(q, q', q+,q-) — M(q, q-,q + ,q'), where 



M(q,q',q+,q-) = -ia J d 4 x J d 4 y^ ql {x)-f^ q (x) 

xD»»{x-y)* q _(y) 7 „* g+ (y), (3) 

and the exchange term M.{q,q-,q+,q') correspond to 
the left and right diagrams in Fig.[TJ respectively. Here, 
^qi, ^ q+ , ^ q _ denote the laser-dressed lepton states 
and D^{x -y) = J ffi )4 e ik '< x ~v) is the free pho- 
ton propagator. Relativistic units with h = c = e.Q = I 
are used. In the presence of the field, the particles are 
characterized by their dressed four-momentum q^ , which 
is related to the free four-momentum outside the field 
by q^ = j/ 1 + " fc p k^. Here, m is the electron mass, k^ 
the laser wave four- vector, and £ = ^A, with the root- 
mean-square value A of the laser vector potential. 

After Fourier series expansions of the periodic parts in 
the laser-dressed leptonic transition currents in Eq. ((3|) , 
the space-time integrations can be performed. The am- 
plitude adopts the form 



M(q,q',q + ,q.)= £ £ 



S^(q + Nk-q' -q+-q_) 
{q-q' + nk) 2 



N>N n=-oo 

xM»(q,q'\n)M^q + ,q_\N-n), (4) 

where A is the total number of photons absorbed, n is 
the number of photons absorbed at the first vertex (see 
Fig. H]), and M^(q, q'\n), M M (q + ,g_|A — n) are complex 
functions of the particle momenta and laser parameters. 



The threshold photon number Nq following from four- 
momentum conservation is the smallest integer A with 



Nuj' > 4ra* 



(5) 



where u/ is the laser photon energy in the average rest 
frame of the projectile and ro* the dressed electron mass 
defined as m 2 = q 2 = m 2 (l + £ 2 ). The intensity depen- 
dence of the dressed mass influences the threshold; for ex- 
ample, Ao = 6 for the SLAC parameters p° = 46.6 GcV, 
uj = 2.35 eV and £ s» 0.3 [3J, whereas Ao = 5 photons 
would suffice in a weaker field (£ < 0.22). 
The total rate is obtained as 



d 3 q + f d 3 q. 



Tj (2tt)3; (2tt)3 



d 3 q' 1 ^ 
(2^4 V 1/4 

spins 



(6) 



with the interaction time T and a statistical factor 1/4 
due to initial spin averaging and the two identical final- 
state electrons. The spin sum can be converted in the 
usual way into voluminous trace products. The multi- 
dimensional integration in Eq. (|6]) is carried out numeri- 
cally by an appropriate Monte Carlo routine. 

Note that in the weak-field, one-photon limit (( C 1, 
A = 1) our approach reproduces the well-known cross 
section for pair creation by a single 7-photon on an elec- 
tron [22)]. The limit lo — > corresponds to trident pah- 
production in constant crossed fields 23]. As a general 
rule, the rate for an A-photon process in a laser field 
scales as £ 2N in the perturbative regime (£ <C 1). 

We emphasize that our laser-dressed QED approach to 
multiphoton trident pair creation incorporates the direct 
process as well as the two-step mechanism [see Eqs. ((T|), 
(J2J], Due to the presence of the laser field, the inter- 
mediate photon may reach the mass shell and become 
real (kl 2 — 0), provided A > 2 holds. The second-order 
Feynman graph then decomposes into two first-order di- 
agrams which precisely correspond to the multiphoton 
Compton and Breit- Wheeler sub-processes 

For on-shell photon momenta k' 2 = (q — q' + nk) 2 = 0, 
as seen in Eq. ((4]) , the photon propagator formally di- 
verges which is mainly due to the infinite spatio-temporal 
extent of the laser field assumed. A regularization proce- 
dure transforms the divergences into finite physical res- 
onances. These resonances have been studied most ex- 
tensively with respect to laser-assisted M0ller scattering 
5]. In our case we regularize the propagator by taking 
the finite interaction time T into account, according to 
l/k' 2 -> l/(fc' 2 + is) with e = 2\k'°\/T. In the SLAC 
experiment [3J, the electron passes through the laser fo- 
cus in roughly T rj 40 fs (focal waist size ~ 4/xm, beam 
crossing angle 17°). This (lab-frame) interaction time 
will be used in the following [25| . 

Fig. [2]shows the trident pair production rate as a func- 
tion of the laser intensity parameter £. For the electron 
momentum and laser frequency chosen, the £ dependence 
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FIG. 2: (color online) Positron rate dependence on £ in the 
head-on collision of a 46.6 GeV electron with a linearly polar- 
ized 527 nm laser wave. The blue dots denote our numerical 
results from Eq. ((SJ. The green short-dashed and red solid 
lines, respectively, are £ 10 and £ 8 power-law fits to the data, 
and the black long-dashed line is an exponential fit of the 
form exp[— 6n/£\ with 6 = 1.73. No = 6, ... ,17 indicate the 
minimum numbers of laser photons in different £ regimes [see 
Eq. J5J]. The inset shows an enlargement of the nonperturba- 
tive f ~ 1 domain on a linear scale. 



gradually changes from a £ 10 behavior (£ < 0.3) to a flat- 
ter £ 8 increase (£ ~ 0.5), and eventually leads into an 
exponential dependence (£ ~ 1) similar to the famous 
Schwinger rate |4] and marking the transition into the 
fully nonperturbative regime. Here, photon orders up to 
N ~ 50 give significant contributions to the total rate. 

The SLAC experiment found a rate scaling of R ~ £ 10 
around £ w 0.3 |3], appearing indicative of the typical 
£ 2JV ° dependence in the perturbative domain. However, 
as mentioned above, Nq = 6 in the SLAC case. In 
fact, our simulation for £ = 0.3 reveals that, on aver- 
age, 6.44 photons are absorbed in total, with 1.62 (4.82) 
photons being absorbed at the first (second) vertex in 
Fig.[T] Hence, in contrast to the common interpretation, 
the SLAC experiment did not operate in the perturbative 
£ 2JV ° -domain but rather observed the onset of nonpertur- 
bative effects which would become more pronounced for 
£ ks 1 (see Fig. [2J. The continuation of the £ 2Ar ° behavior 
at £ <C 1 into the intensity domain where this Np-th or- 
der channel has already closed was discussed in [26j and 
shown for the example of strong-field atomic ionization. 

From Fig. [5] a total lab- frame positron rate per projec- 
tile of « 4 x 10 4 s _1 results by averaging over a Gaus- 
sian laser focus peaked at £ = 0.3. It is in reasonable 
agreement with the experimental result, where about 100 
e + were created in 22,000 shots from 10 7 electrons in the 
interaction region [3j. The corresponding rate is slightly 
smaller, -R cxp ~ 10 4 s _1 , which is mainly attributable to 
shot-to-shot intensity variations (0.2 < £ pea k ^ 0.3). 

The energy distribution of the positrons is displayed 
in Fig. [3] Our calculations for the SLAC parameters [3| 
reproduce the measurements very well (apart from the 
data point at ~ 15 GeV), see Fig.[3Ja). In addition we 
find that the spectral maximum remains at « 13 GeV 
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FIG. 3: (color online) Positron energy spectra in electron 
collisions with a 527 nm laser beam. (a) £ = 0.3 and 
p° — 46.6 GeV; the dots with an error bar are the experi- 
mental data from [jj. The three calculated curves refer to 
different center-of-mass energies, corresponding to different 
beam crossing angles as indicated, and have been normalized 
to the same height to facilitate their comparison, (b) Non- 
perturbative domain at £ = 1 and p° — 17.5 GeV; the blue 
long-dashed (red short-dashed) line refers to the partial spec- 
trum including up to 30 (35) absorbed photons. The black 
solid line shows the full spectrum (N < 45). 



when the collision angle is varied from 0° to 28°. This 
explains why the peak position in the experiment is not 
blurred by averaging over the laser focus and electron 
beam profile. We note besides, that the mean energy of 
the recoiled electron after head-on collision is 16.8 GeV, 
corresponding to an energy loss of 64%. 

The fully nonperturbative regime could be probed, 
e.g., by utilizing 17.5 GeV electrons from the upcom- 
ing European XFEL beamline at DESY (Hamburg, Ger- 
many) [27J combined with a table-top 10 TW laser sys- 
tem. Here, the contributions from many photon orders 
form the e + spectrum, see Fig.[3]Jb). Such an experiment 
would represent a non-standard application of the XFEL 
electron beam, usually serving to generate x-ray light. 

In the SLAC experiment, the pair production was 
found to be largely dominated by the two-step mech- 
anism |3j, leading to a process rate proportional to T 
[H HH. We note that the direct process (TJ could be 
revealed as well by another suitable choice of collision 
parameters. In the laboratory frame, the threshold for 
the direct process with A-photon absorption is found to 
be wbh ^ jj — , which lies below the corresponding two- 
step mechanism threshold cjbw ^ 2{V7j~i) W ^ ^ 
being the projectile Lorentz factor. An example is shown 
in Fig. [4] where a 17.5 GeV electron 27J collides with 
an intense soft VUV pulse [2jj]. Below wbw ~ 18 eV, 
two-photon pair creation is possible via the direct mech- 
anism only, which can be measured separately at VUV 
intensities ~ 10 13 W/cm 2 (£ = 10~ 4 ). In this region, 
the two-step mechanism requiring an additional photon 
is strongly suppressed. 

The SLAC experiment relied on the high-energy elec- 
tron beam from a large-scale linear accelerator. Nowa- 
days, corresponding pair creation studies could be per- 
formed with compact laser wakefield accelerators pro- 
ducing few-GeV electron beams 30]. Assuming a laser- 
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FIG. 4: Laser frequency dependence of the pair creation rate 
in the head-on collision of a VUV pulse (£ = 10~ 4 ) with a 
17.5 GeV electron. Shown are the separate contributions from 
the direct (circles) and two-step (triangles) processes whose 
sum yields the total rate. 
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FIG. 5: Relation of the optical (527 nm) laser intensity 
and the electron energy to give an observable positron rate 
~ 10 5 s _1 in the lab frame. The pair creation mechanism 
is changing from the perturbative few-photon regime to the 
quasi-static (i.e. low-frequency, high-intensity) regime, as in- 
dicated. 

accelerated 5 (2) GeV electron colliding with a second op- 
tical laser pulse of intensity ~ 10 20 (10 21 ) W/cm 2 , an ob- 
servable pair creation rate of ~ 10 5 s _1 in the manifestly 
nonperturbative regime £ ~ 3 (8) results [3l| . At the en- 
visaged high-power (~ 10 25 W/cm 2 , £ ~ 10 3 ) facility ELI 
[32j ]. comparable rates can be expected for p° ~ lOMeV 
already. The pair creation occurs in a quasi-static regime 
[33l | here since, in the electron frame, the low-frequency 
laser field resembles a constant crossed field 23[ on the 
time scale necessary for the field to create a pair. Another 
all-optical scheme for pair creation in two laser beams 
emplo ys a seed electron being accelerated directly by the 
fields [H . 

Fig. [5] provides an overview of our results on the total 
pair production rates in various interaction regimes. 

In conclusion, a complete QED treatment of multi- 
photon trident pair creation in electron-laser collisions 
has been presented. It was shown that the SLAC exper- 
iment 0] observed nonperturbative signatures. Future 
experimental studies, allowing for a disentanglement of 
the competing production processes in the perturbative 
regime and probing the transition to the fully nonper- 
turbative domain, may rely on relativistic electrons from 
XFEL beamlines or laser accelerators. 
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